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ABSTRACT 

We present novel evidence for a fine structure observed in the net-circular polariza- 
tion (NCP) of a sunspot penumbra based on spectropolarimetric measurements utiliz- 
ing the Zeeman sensitive Fe I 630.2 nm line. For the first time we detect a filamentary 
organized fine structure of the NCP on spatial scales that are similar to the inhomo- 
geneities found in the penumbral flow field. We also observe an additional property 
of the visible NCP, a zero-crossing of the NCP in the outer parts of the center-side 
penumbra, which has not been recognized before. In order to interprete the observa- 
tions we solve the radiative transfer equations for polarized light in a model penumbra 
with embedded magnetic flux tubes. We demonstrate that the observed zero-crossing 
of the NCP can be explained by an increased magnetic field strength inside magnetic 
flux tubes in the outer penumbra combined with a decreased magnetic field strength 
in the background field. Our results strongly support the concept of the uncombed 
penumbra. 

Subject headings: Sun: sunspots — magnetic fields — photosphere 



1. Introduction 

The asymmetry of the Stokes parameters with wavelength is a powerful diagnostic tool to 
probe the relation between flows and the magnetic field in the atmosphere of sunspots. Line 
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asymmetries convey information about lateral (within one resolution element) and the line-of- 
sight (LOS) gradients and discontinuities of physical parameters. Since the first systematic char- 
acterizatio n and interpretation of these asy mmetries and their spatial distribution throughout a 
sunspot by ISanchez Almeida & LitesI (119921 and references therein) the complementary improve- 
ments achieved in both observational techniques and theoretical modeling, have advanced our 
understanding of particularly the penumbral inhomogeneities. It is general consent that a combi- 
nation of non-trivial magnetic field and flow topologies along the LOS are the cause of the ob- 
served penumbral Stokes asymmetries. A simple and instructive way to characterize the asymme- 
try of Stokes V profiles is through the net-circular polarization (NCP), defined as ,jY = JV{X)dX, 
where the integral is performed over one spectral line only. The sum of the NCP of individ- 
ual lines than lead s to the well-known broad-band circular polarization (BBC?) first noticed by 
lUing et al.l (|l974al]br) . The appearance of the NCP is very distinct in the penumbra and shows a 
peculiar difference when between observations in the visible or in the near infrared wavelength 
region: in the visible at Fe 1 630. 2 nm the NCP maps tend to be symmetric about the disk cen 
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ter and sun center connection line ([Sanchez Almeida & Litesll 19921: IWestendorp Plaza et al.ll2001 



2002h while infrared NCP maps at Fe 1 1564.8 nm reveal an antisymmetr i c dist ribu 



tion Jschlichenmaier & Colladosll2002l : lBellot Rubio et al.ll2004h . Ischhchenmaier et al.l(l2002h and 
Miiller et ah! ( 2002 ) demonstrate that the un combed penumbral atmosph ere ( Solanki & Montavon 
1993|) mimicked by the moving tube model dSchlichenmaier et al.lll998h reproduces the azimuthal 
behavior of the NCP and, furthermore, convincingly show that the incisive discrepancy between 
the visible and the near infrared can b e understood when the effects of anomalous dispersion 
([Landolfi & Landi deglTnnocentilll996|) are included. 



So far penumbral NCP map s of the Fe 1 630.2 nm line have been only presented before b y 
Sanchez Almeida & LitesI (119921) (in form o f a contour line map), IWestendorp Plaza et al.l (|200ll) . 
Miiller et al.l (|2002r) and iMuller et al.l (|2006l) . These maps visualize very well the symmetric be- 
haviour of the NCP but lack the spatial resolution to reveal the presence of any potential fine 
structure. 

In this Letter we focus on high-spatial resolution spectropolarimetric measurements of a 
sunspot to investigate inhomogeneities in the penumbral magnetic field and flow field. Instead 
of utilizing inversion techniques we calculate the NCP and line-of-sight (LOS) velocities and elab- 
orate on the properties of the same. In order to interpr et the observational results we use the 3D 
geometric flux tube model VTUBE (IMuller et al.ll2006l) and perform subsequent radiative transfer 
calculations from which NCP and Doppler maps can be synthesized. 
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2. Observations and Data Reduction 



The spectropolarimetric observ ations were obtained on 2005 Ju ne 30, with the Diffraction 
Limited Spectro-Polarimeter (DLSP, Sankarasubramanian et al. 20041) [^located at the NSO Dunn 
Solar Telescope (DST), Sacramento Peak, New Mexico. The DLSP was operated in its high- 
resolution mo de (0.^0893 pixe l ^ along the slit) and in conjunction with the high-order adaptive 
optics system ( Rimmelj2004 ). The slit width of 12 microns is matched to the detector image scale 
along the spectrograph slit. In spectral direction the instrument samples the solar spectrum from 
630.0nm to 630.4 nm with 2.1 pmpixel^^. A two-dimensional map was generated by moving 
the slit across the leading sunspot (positive polarity) of AR NOAA 10781 located at jU = 0.74 
(corresponding to a viewing angle of 42°) with a step size of Of'089 covering a total field-of- 
view of 58^'9 (660 steps). The slit was oriented along solar North-South and the scan direction 
of the DLSP was solar West-East. At each slit position 16 images of 30 ms exposure time each 
were accumulated for each of the four individual modulation states leading to a total integration 
time of 2 sec. A full map is accomplished within 56min. The measurements are compensated for 
polarization e ffects caused by the teles cope and the polarimeter following the calibration procedure 
developed by ISkumanich et al.l (1 1997b . The very good seeing in conjunction with the high-order 
AO system led to exceptional stable observing conditions while scanning the solar image. We 
estimate our effective spatial resolution in the DLSP map to be about 0''4. 

We calculate maps of the net circular polarization by integrating the Stokes V signal over 
the Fe 1 630.25 nm line. Because of its relatively minor contamination by one of the two telluric 
lines the Fe 1630.15 nm line is used to deriv e line-of-sight (LOS) Doppler velocities. We apply 
a Fourier phase method (ISchmidt et al.ll 19991) to the I + V and I — V profile of the Fe 1 630. 1 5 nm 
line to calculate line shifts. An absolute velocity calibratio n is performed by taking advantage of 
the telluric O2 lines (see e.g. iMartinez Pillet et al.lll997h present in the recorded spectrum. In 
contrast to the astrophysical convention we assign a negative sign to redshifts and a positive sign to 
blueshifts: dark areas in the Dopplergram move away from the observer, while bright areas move 
towards the observer. 



3. Results 

Well-reported properties of the penumbral NCP (see Sanchez Almeida & Litesll992 : Solanki & Montavon 



1993|) are the sign change along azimuthal cuts centered around the spot center and an increased 
magnitude on the limb side. Based on spectropolarimetric data with much higher spatial resolution 



^The DLSP was built in collaboration with the High- Altitude Observatory. 
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than utilized before (Figure [T]) we additionally find (a) evidence for a fine-structure in the NCP 
and (b) a sign-reversal of the NCP as a function of radial distance from the spot center on the 
center- side penumbra which has not been recognized before. 

Confirming previous results we find that the NCP is symmetric in the inner and middle penum- 
bra w.r.t. the line that connects the disk center with the spot center. The same property applies to 
the LOS velocity map which indicates an intricate coupling between penumbral flows and the sign 
of the NCP. This is corroborated by the fine structure we detect in the NCP. 

The spatial distribution of the NCP reveals filamentary organized inhomogeneities through- 
out the penumbra with very localized enhancements of elongated shape on a spatial scale of 0'.'3. 
Although we detect these enhancements with no preferred location inside the penumbra those with 
the largest value of the NCP appear on the limb-side penumbra in proximity to where the magne- 
togram (not shown) indicat es an apparent neutral line. They are related to anomalous (multi-lobe) 



Stokes V profiles (see also ISanchez Almeida & Liteslll992h and thus do not necessarily show up 



in the flow signal (Figure [H right). Anomalously shaped Stokes V profiles indicate either spa- 
tially unresolved contributions to the resolution element or LOS effects of at least two differently 
inclined magnetic components with opposite polarity that are shifted w.r.t. each other and viewed 
upon under an oblique angle. The NCP on the center-side penumbra is correlated with normal 
2-lobe antisymmetric Stokes V profiles. On smaller spatial scales related to the flow channels 
(FigurefH right) that carry most of the Evershed flow we are unable to state conclusively that there 
is a good overall correlation between the NCP and the flow signatures. Figure [2] exemplifies this 
statement along three azimuthal cuts through the penumbra representative for three different radial 
distances from spot center (as indicated by the ellipses in Figure [H left). 

Furthermore, we detect for the first time a difference in the symmetry properties of the NCP 
observed between the outer limb- and center-side penumbra, namely a sign reversal of the NCP in 
the outer center-side penumbra. We verify the visual impression of the NCP zero-crossing on the 
center-side penumbra in Figure [3] by means of the radial dependence of the NCP when averaged 
over the limb-side (90°-270°) and center-side penumbra (270°-90°) separately. The limb-side 
segment without penumbra does not contribute to the averaging. 

A comparison of the two curves clearly demonstrates the differences between the NCP on the 
limb- and center-side penumbra in radial direction. On the limb side the NCP has one sign only 
throughout the penumbra while on the center side the NCP reverses its sign at a radial distance 
of ~0.757?, (where R denotes the sunspots radius) peaks at 0.95i? and then falls off outside the 
visible penumbra. It should be noted that the course of the radial variation of the NCP and the 



^Beyond the penumbra crossing the visible boundary of the sunspot the NCP shows a motley behavior which we 
attribute to the presence of moving magnetic features that are known to show asymmetric Stokes profiles. 
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exact location of the zero-crossing point within the penumbra depends also on the calibration for 
instrumental polarization, i.e. predominantly the crosstalk from Stokes / to V. We estimate that our 
calibration for this particular crosstalk is better than 0.1 %. The NCP averaged over areas where 
we see no significant Stokes V signal is in the range of 1.5x lO^^pm. Significant changes of the 
radial variation of the NCP, however, take place only if the residual crosstalk from / to V is in the 
range of 1 % and larger. 



4. Model Calculations 

In order to understand the observed sign reversal of the NCP, we have carried out radiative 
transfer calculations from which synthetic maps of observable quantities are obtained. At first 
sight, one may conjecture that the observed radial sign reversal of the NCP is caused by downflows 
in the outer penumbra. While this is indeed a possible solution that is compatible with the observed 
NCP maps, the Doppler maps constrain the details of the model atmosphere even further. 

It turns out that a model atmosphere with strongly downwards inclined flux tubes in the outer 
penumbra is not able to explain simultaneously both the observed sign reversal of the NCP and 
the constant sign of the Doppler maps on the spot's center side. However, we find that both fea- 
tures can be reproduced by magnetic flux tubes that have a magnetic field strength that is stronger 
than the background magnetic field in the outer penumbra while it is weaker in the inner part as 
visualized in Figured (top left). The assumption is not unrealistic as two-component inversions 
of s unspots indicate a higher field strength in the tube component in the extreme outer penumbra 



(e.g. iBellot Rubio et al.ll2004l) . 



Figure |4] shows the synthetic NCP (middle) and Doppler (right) maps for a penumbral at 



mosphere, calculated with the VTUBE model of iMuller et al.l (|2006h . We adopt a very generic 



atmospheric configuration with horizontal magnetic flux tubes of 50 km ra dius with a flow ve 



locity of lOkm/s, embedded in a radially symmetric background field as in iMuller et al.l (|2006l) . 
We assume a magnetic field strength inside the flux tubes that increases linearly from the inner 
penumbra (where it is 500 G weaker than the background field) to the outer penumbra (where it is 
500 G stronger than the background field). It turns out that this configuration is quite robust, and 
the details of the magnetic field strength affect primarily the radial location of the zero-crossing 
of the NCP. The axes of the flux tubes are placed at 30 km above T500 = 1 • Raising the flux tubes 
mostly increases the absolute value of the NCP and would require a lower filling factor to explain 
the observed NCP values. 



The effect of a jump in the magnetic fi eld strength along the LOS on the NCP h as been studied 
for a simpler atmospheric configuration by lLandolfi & Landi deglTnnocentil(|l996|) who derived a 
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sign rule for the NCP due to this so-called AS effect. This expression is proportional to the differ- 
ence in the magnetic field strengths between the static background and the part of the atmosphere 
containing the flow. As a consequence the NCP changes sign when AS changes sign. When com- 
paring this rule with both observational data and more realistic models, one should note that the 
Afl effect is but one of several NCP generating effects at work. Therefore, the locations of the sign 
reversals of AB and the NCP do not necessarily coincide. In particular the sum of the AB effect 
and changes in the inclination (Ay effect) results in the different behaviour of the NCP on the limb 
and center side: on the limbside the NCP decreases with radial distance because of the increasing 
AB effect but does not change the sign because the Ay effect con tributes positive to the NCP. One 



should also bear in mind the assumptions on w hich the findings oflLandolfi & Landi deglTnnocenti 



(Il996|) are based. A discussion of the results o flLandolfi & 



text of flux tubes in the penumbra is given by lMuUer 



andi degrinnocentil (1 1996h in the con 



and 



MulleretaL 



mm. 



Our model solution is in qualitative agreement with both the observed NCP and the observed 
LOS velocity. It should be pointed out though that while the simultaneous agreement between the 
NCP and Doppler maps places a strong constraint on the possible atmospheric models our solution 
may not necessarily be unique. 



5. Discussion and Conclusions 



In this paper, we have presented novel evidence for a spatial fine structure in the visible NCP 
of the Fe 1 630.25 nm line. We find for the first time that the NCP is structured in radial filaments 
very mu ch alike maps of other line parameter s, e.g. the intensity, equi valent width and line width 
(see e.g. Ijohannessoni 1 19931 : iRimmeld 1 1995alJbl : iTritschler et al.l 120041) . The spatial scales of the 
NCP filaments are similar to those found in the observed LOS velocity map although we cannot 
find a good overall correlation between the signals of the NCP and the LOS veloc ity. We conjecture 
that higher spatial resolution is needed to settle this issue. iMartinez Pilleti (|1997h find a correlation 
between the azimuthal variation of the NCP and the magnetic field inclination on angular scales of 
10°. Although we cannot provide information about the inclination at this point we find that our 
observations show azimuthal variations of the NCP that are only of the order of 4°. 

Furthermore, we have observed a zero-crossing of the NCP on the outer center-side penumbra 
and have demonstrated that this signature can be explained by an increased magnetic field strength 
inside magnetic flux tubes in the outer penumbra, while the magnetic field strength is decreased 
in the background field (A5-effect). To this end, we have compared high-resolution observations 
with maps of the NCP and the LOS velocity that have been synthesized , from a 3D georn etric 
sunspot model that has been coupled with a ID radiative transfer code (see lMuUer et al.ll2006l for 
details). We point out that the space of possible atmospheric models that can reproduce the ob- 
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servations is significantly constrained by taking into account the info rmation contained in Dopple r 
maps in addition to th e NCP maps. This p aper extends the work by ISolanki & MontavonI (| 19931) . 



Martinez Pilletl (|2000|) . iMuUer et all (|2002h . and lMuUer et all (|2006h and clearly demonstrates the 



need to model the "uncombed" penumbra of sunspots in greater detail. This would in clude also tak- 



ing int o account the effect of field line curvature around flux tubes as considered by iBorrero et al. 
J2OO7I) . 



REFERENCES 

Bellot Rubio, L. R., Balthasar, H., & CoUados, M. 2004, A&A, 427, 319 
Borrero, J. M., Bellot Rubio, L. R., & Mueller, D. A. N. 2007, ArXiv e-prints, 707 
lUing, R. M. E., Landman, D. A., & Mickey, D. L. 1974a, A&A, 35, 327 
— . 1974b, A&A, 37,97 
Johannesson, A. 1993, A&A, 273, 633 

Landolfi, M., & Landi deglTnnocenti, E. 1996, Sol. Phys., 164, 191 

Martinez Pillet, V. 1997, in ASP Conf. Ser. 118, Advances in Physics of Sunspots, ed. 
B. Schmieder, J. C. del Toro Iniesta, & M. Vazquez, 212 

Martinez Pillet, V. 2000, A&A, 361, 734 

Martinez Pillet, V., Lites, B. W., & Skumanich, A. 1997, ApJ, 474, 810 

Miiller, D. 2001, Polarisation von Linien im Spektrum einer Sonnenfleckenpenumbra. Diploma 
Thesis, University of Freiburg 

Miiller, D. A. N., Schlichenmaier, R., Fritz, G., & Beck, C. 2006, A&A, 460, 925 

MuUer, D. A. N., Schlichenmaier, R., Steiner, O., & Stix, M. 2002, A&A, 393, 305 

Rimmele, T. R. 1995a, A&A, 298, 260 

— . 1995b, ApJ, 445,511 

— . 2004, ApJ, 604, 906 

Sanchez Almeida, J., & Lites, B. W. 1992, ApJ, 398, 359 



-8- 

Sankarasubramanian, K., GuUixson, C, Hegwer, S., Rimmele, T. R., Gregory, S., Spence, T., et 
al. 2004, Proceedings SPIE, 5171, 207-218 

Schlichenmaier, R., & CoUados, M. 2002, A&A, 381, 668 

Schlichenmaier, R., Jahn, K., & Schmidt, H. U. 1998, ApJ, 493, L121 

Schlichenmaier, R., Miiller, D. A. N., Steiner, O., & Stix, M. 2002, A&A, 381, L77 

Schmidt, W., Stix, M., & Wohl, H. 1999, A&A, 346, 633 

Skumanich, A., Lites, B. W., Martinez Fillet, V., & Seagraves, R 1997, ApJS, 1 10, 357 
Solanki, S. K., & Montavon, C. A. R 1993, A&A, 275, 283 

Tritschler, A., Schlichenmaier, R., Bellot Rubio, L. R., the KAOS Team, Berkefeld, T., & Schelenz, 
T. 2004, A&A, 415, 717 

Westendorp Plaza, C, del Toro Iniesta, J. C, Ruiz Cobo, B., Martinez Pillet, V., Lites, B. W., & 
Skumanich, A. 2001, ApJ, 547, 1130 



This preprint was prepared with the AAS liT^ macros v5.2. 



-9- 




42° 



- w 




Fig. 1. — Continuum intensity (left), net circular polarization (middle) and LOS velocity (right) 
of AR NOAA 10781 observed on June 30, 2005, between 14:24 UT and 15:03 UT with the DLSP. 
The Dopplergram is calculated from the Fe 1 630.1 5 nm line while the NCP is determined from 
the Fe 1 630.25 nm line. The arrow points towards disk center (DC). The three ellipses indicate 
azimuthal cuts corresponding to radial distances from the center of the umbra of 0.5R, O.IR, and 
0.97? in units of the spot radius, R. The contour lines mark the inner and outer boundary of the 
penumbra as determined from thresholding a smoothed version of the continuum map at 0.70/c 
and 0.96/c in units of the continuum intensity of the quiet sun. 
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Fig. 2. — Azimuthal variation of the NCP (black) and LOS velocity (grey) measured counterclock- 
wise from disk center direction for three different radial distances from spot center as indicated in 
Figure [T] (left). The NCP is calculated from the Fe 1 630.25 nm line while the LOS velocity is 
derived from the Fe 1630.15 nm line. The shaded areas where no data points are plotted (141.5°- 
225.0°) mark the angle segment where the penumbra is missing. Horizontal dotted lines indicate 
zero NCP and LOS velocity. 




Fig. 3. — Radial dependence of the mean LOS velocity (< vlos >) and the NCP averaged over 
the limb-side (< . >limb) and center-side (< . >center)5 respectively. The penumbra missing 
segment (141.5°-225.0°) is not contributing in the averaging process. The vertical dotted lines 
mark the inner and outer penumbral boundary. Horizontal dotted lines indicate zero NCP and 
LOS velocity. The radial distance is given with respect to the spot radius. 
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Fig. 4. — Left : Radial variation of the magnetic field strength inside the flux tube and in the back- 
ground field (top) and variation of the NCP on the center- and limb-side penumbra (bottom). Syn- 
thetic maps of the net circular polarization (middle) and LOS velocity (right) of the Fe 1 630.25 nm 
line for a sunspot at 6 = 42°. To illustrate the filamentary structure of the penumbra, we alternate 
between flux tubes with a radial width of 2° and "gaps" of stationary background atmosphere of 
the same width. 



